Metabolic Modeling of T-cells
Sara Cardoso
1 From gene expression to reaction presence
When looking into the UMAPs for scRNAseq data using all genes versus using only the metabolic genes, there is a marked difference in the capacity to separate the cells into their respective cell-types. In fact, only the proliferative T-cells are clearly separated from the rest when using only metabolic genes, which are very well known for their drastic change in metabolism to adjust to their proliferative demands.
Regarding the pseudo-bulk data, which was used to create a cell-type (if present) model for each sample in each individual, we can still see the same separation of the proliferative T-cells. Furthermore, other cell-types seem to be grouping together.
Next, machine learning with a random forest classifier was performed using either pseudo-bulk data or reaction presence data (which represents the reactions that were considered present when constructing the cell-type specific models).
The evaluation metric used to assess the prediction capability was the Mathews correlation coefficient (MCC), as it is appropriate to classify multi-class problems, perfectly symmetric (no class is more important than the other), and not sensitive to class-imbalance (i.e., when the different classes are not evenly represented, which happens in this data, with proliferative T-cells having far less models than other cell-types). The general MCC formula is the following:
\[ MCC = \frac{TP * TN - FP * FN}{\sqrt{(TP + FP)(TP + FN)(TN + FP)(TN + FN)}} \]
The MCC metric varies from -1 to 1. An MCC value 1 means that all cell-types were correctly classified, while a value of -1 means that all cell-types were not well classified. An MCC value of 0 means that the classifier is no better than random guessing.
The dataset that leads to the best results is the pseudo-bulk data, with an MCC greater than 0.75. Prediction capacity slightly decreases when using reaction presence data, with an MCC of around 0.56. This shows that the models do a good job at representing the cell-types at the level of what reactions should be included or not to accurately represent the transcriptomics data.
2 Pathways’ presence
We next checked the percentage of reactions present in each metabolic pathway (relative to all reactions of that pathway in the generic human model) for all models of each cell-type, and visualised the results in an heatmap.
Several metabolic pathways are highly present across all cell-types, namely: aminoacyl-tRNA biosynthesis, arachidonic acid metabolism, estrogen metabolism, fatty acid activation (cytosolic), keratan sulfate biosynthesis, leukotriene metabolism, linoleate metabolism, protein degradation, and vitamin E metabolism.
There are other metabolic pathways that are present somewhat to an equal extent across all cell-types: butanoate metabolism, cholesterol biosynthesis 1 (Bloch pathway), chondroitin / heparan sulfate biosynthesis, glycerolipid metabolism, histidine metabolism, lysine metabolism, porphyrin metabolism, O-glycan metabolism, oxidative phosphorilation, serotonin and melatonin biosynthesis, vitamin D metabolism, and xenobiotics metabolism.
For the remaining pathways, some differences between cell-types can be observed. In fact, overall, memory and regulatory T-cell models have more pthway presence than other cell-types.
Biotin metabolism, for example, has a small reaction presence in IL17+ CD4 T-cells, as well as in proliferative CD4 T-cells, when compared to naive, regulatory, follicular and memory CD4 T-cells. In fact, biotin deficiency was shown to enhance secretion of pro-inflammatory cytokines from CD4 T-cells expressing IL17, as well as differentiation of naive and memory CD4 T-cells toward \(T_{H1}\) and \(T_{H17}\) cells. On the other hand, this deficiency decreases differentiation toward anti-inflammatory regulatory T-cells (Elahi et al. 2018). Finally, biotin deficiency was found to lead to a decrease in CD4 T-cell proliferation (Elahi et al. 2018). In fact, if we plot the presence of biotin metabolism pathway in relation to biomass production in proliferative CD4 T-cell models, we can see a slight positive correlation:
For regulatory CD4 T-cell models, we can see that the reaction presence in some pathways is way lower in models from normal matched tissue than in those from their tumour counterparts. The following heatmap shows this:
In this heatmap, we can see that biotin metabolism is clearly more present in regulatory CD4 T-cell models from tumour tissues. Considering that biotin deficiency decreases differentiation toward anti-inflammatory regulatory T-cells (Elahi et al. 2018), it makes sense that this pathway would be more present in regulatory CD4 T-cells from tumour tissue, where they might be exerting a highly immune suppressive role that helps tumour growth, while in the normal matched, non-inflammatory, tissue it is not as necessary.
The other pathways that have a clear separation between tumour and normal matched mucosa derived models are Biopterin metabolism, Bile acid synthesis, and Fatty acid biosynthesis.
We can see something similar for naive CD8 T-cells, with pathways Glycosylphosphatidylinositol (GPI)-anchor biosynthesis, keratan sulfate degradation, and Beta-oxidation of fatty acids:
However, this time, the metabolic pathways are more present in naive CD8 T-cell models derived from normal matched mucosa tissues than those from tumour tisues.
Fatty acid related pathways are also worth looking into:
Most models from memory CD4 and CD8 T-cells have a high reaction presence of \(\beta\)-oxidation and biosynthesis of FAs, as seen in literature (O’Sullivan et al. (2014)).
Regarding regulatory CD4 T-cells, most models have a high reactions presence of \(\beta\)-oxidation of FAs, also corroborated by the literature (Cluxton, Moran, and Fletcher (2019)). However, regulatory T-cells have been found to not entirely rely in FA oxidation (Cluxton, Moran, and Fletcher (2019), Raud et al. (2018)). This can be the case with those models, all from normal tissue samples, that present low reaction presence of \(\beta\)-oxidation of unsaturated FAs, phytanic acid and even-chain FAs in the peroxisome. These same models also have low reaction presence of odd-chain FA biosynthesis.
Nevertheless, all models from memory CD4 and CD8 and regulatory CD4 T-cells have high reaction presence of FA activation in both cytosol and endoplasmic reticulum.
For models of IL17+ CD4 T-cells, those from normal matched tissue have really low reaction presence of all \(\beta\)-oxidation pathways, while having high reaction presence in the biosynthesis pathway of unsaturated FAs, which is corroborated in the literature (Berod et al. (2014)). On the other hand, models from tumour tissue have high reaction presence of some of the \(\beta\)-oxidation pathways, more specifically oxidation of poly-unsaturated FAs in mitochondria and unsaturated FAs (n-7) in peroxisome.
Regarding naive T-cells, they are characterised by high FA oxidation (MacIver, Michalek, and Rathmell (2013)). In fact, we can see this across most naive CD4 T-cell models, while for naive CD8 T-cell models only those from normal matched tissue show high presence of FA oxidation.
In general, reaction presence of pathways related to fatty acids are in agreement with the literature, with the exception of IL17+ CD4 and naive CD8 T-cell models from tissue samples. This could mean that the tumour micro-environment does affect IL17+ CD4 and naive CD8 T-cells’ FA related metabolism.
3 Biomass vs ATP production
As expected, proliferative CD4 and CD8 T-cells seem to have a higher biomass flux than their naive counterparts. Regarding the remaining cell-types, each seems to have a varied biomass flux across the different models, with the exception of cytotoxic CD8 and IL17+ CD4 T-cells, which have relatively smaller biomass fluxes.
Naive CD8 T-cells have markedly more ATP production than their proliferative counterparts. However, the same thing for naive vs proliferative CD4 T-cells is not clear.
When comparing biomass and ATP production, ATP production is clearly bigger in IL17+ CD4 T-cells and naive CD8 T-cells. In proliferative CD4 and CD8 T-cells, the flux going through the biomass reaction is clearly bigger than their ATP production.
We then checked how different biomass and ATP production were between cell-type models from normal vs tumour tissue:
Biomass
ATP Production
For cytotoxic CD8 and regulatory CD4 T-cells, models from tumour tissues have higher biomass flux than those from normal matched mucosa. For memory CD8 T-cells, there is increased biomass in models from tumour border.
4 Sources of FADH2 and NADH
In general, all cell-types resort obtain most of their FADH2 and NADH from Fatty Acid Oxidation (FAO), which is expected according to the literature.
FADH2
NADH
5 FAs Uptake
In general, cytotoxic CD8, naïve CD4 and proliferative CD4 and CD8 T-cells are the cell-types with the least overall uptake of fatty acids (FAs). Although on average the other cell-types uptake more FAs, the uptake is still very varied across the models of each cell-type.
This is in line with the literature. Proliferative T-cells rely more in FA synthesis and decrease FA oxidation relative to their naïve counterparts (Dumitru, Kabat, and Maloy (2018), Kolan et al. (2020)), which might result in decreased necessity for uptake of FAs. Memory CD8 T-cells promote FA oxidation, unlike the activated ones (Kolan et al. (2020)), which might explain why FA uptake for memory CD8 T-cells is generally higher than cytotoxic and proliferative CD8 T-cells.
We then checked how different biomass and ATP production were between cell-type models from normal vs tumour tissue:
6 Medium Changes
6.1 No Tryptophan
[B2]
It is expected that absence of tryptophan from medium causes T-cells’ biomass to decrease, as it has been shown that indoleamine-pyrrole 2,3-dioxygenase (IDO), which catalyses tryptophan metabolism in the kynurenine pathway, inhibits T-cell activation by tryptophan deprivation and by promoting the expansion of regulatory T-cells (Le Bourgeois et al. (2018)).
Indeed, the biomass of all models, across all cell-types, decreases to zero or very close to zero once tryptophan is removed from medium.
6.2 No Oxygen
[B4, B16]
In general, it is expected that T-cells suffer reduced proliferation in an environment with no oxygen (Dumitru, Kabat, and Maloy (2018), Loeffler, Juneau, and Masserant (1992)), even though reduced amounts of oxygen up-regulate genes involved in glycolytic ATP production and down-regulates the OXPHOS pathways (Tripmacher et al. (2008)), associated with higher (Carswell, Weiss, and Papoutsakis (2000), Haddad et al. (2004)) or no effect (Tripmacher et al. (2008)) in proliferation. Thus, it could be expected that removing oxygen from the metabolic models’ medium would result in decreased biomass. However, that is not the case in any of the cell-types, even though very few models, across most cell-types, do suffer a decrease in biomass flux. A cell-type where no models suffer a decrease in biomass when oxygen is removed from the medium is IL17+ CD4 T-cells. In fact, hypoxia is associated with inducing Th17 cell proliferation (Krzywinska and Stockmann (2018)).
Still, it has been pointed out that the impact of oxygen in cell viability relies on the type of stimulus that the stimulated cultures received, as two different stimuli revealed different impacts of oxygen levels on T-cells proliferation (Atkuri, Herzenberg, and Herzenberg (2005)). This might help explain why some models show a decrease in biomass while others from the same cell-type show the contrary.
One important thing to take in mind is that metabolic models only capture metabolic pathways, i.e., signaling and/or regulation pathways are not captured in this models. This means that the regulation related to the hypoxic inducible factors (HIFs) are not captured in these models.
Finally, as expected, the reaction catalised by oxygen-dependent prolyl hydroxylase domain enzymes (PHDs) does not have any flux in any model when no oxygen is present in the medium:
Next, we checked if the models that uptaked the most amounts of oxygen were the ones that suffered the most in biomass flux when oxygen was no longer available:
What happens to the biomass of the cell-types from normal tissue vs tumour tissue, when there is no oxygen:
6.3 No Glucose and/or Glutamine
[B8]
Unavailability if glutamine in medium, with glucose present, decreases biomass flux to or very close to zero across all cell-types.
As reported in the literature, glutamine seems to be essential for all T-cells’ proliferation (Dumitru, Kabat, and Maloy (2018), Kolan et al. (2020), Le Bourgeois et al. (2018)), especially because it acts as a nitrogen donor for DNA and RNA nucleotide production (Dumitru, Kabat, and Maloy (2018), Kolan et al. (2020)). Indeed, DNA and RNA production decreases to zero, or very close to zero, when no glutamine is available in the models’ medium.
[B9]
DNA Production
RNA Production
Removing both glucose in addition to glutamine does not seem to make much difference in the biomass flux predicted. In fact, when only removing glucose from medium, biomass only decreases slightly in the cell-types follicular CD4, proliferative CD8, and regulatory CD4 T-cells. In fact, most models in these cell-types do not suffer changes in the biomass flux, with only some suffering from biomass flux decrease. This, however, is not corroborated by literature, as it has been reported decreased T-cell proliferation rates in glucose-deficient media (Dumitru, Kabat, and Maloy (2018)).
For those models whose biomass is affected as expected in the literature, increasing glutamine availability only leads to biomass increase in very few models. In fact, presence of high levels of alternative energy sources like glutamine has been found to not help increasing proliferation in glucose-deficient media (Dumitru, Kabat, and Maloy (2018)).
No Glutamine and Glucose
No Glucose
All together
6.4 No Nucleotides
[B17]
Overall, there is no difference in biomass flux and production of DNA when no nucleotides are available in the medium. This goes in line with a study (Varanasi, Ma, and Kaech (2019)) on in vitro vs in vivo metabolismof CD8+ T-cells, where the effector cells where shown to almost entirely rely on de novo nucleotide biosynthesis.
Biomass
Production of Nucleotides
Production of nucleotides assessed through the DNA production reaction:
